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1. The designated Office is hereby notified of its election made: 

j X | in the demand filed with the International Preliminary Examining Authority on: 

02 June 1999 (02.06.99) 

[ | in a notice effecting later election filed with the International Bureau on: 



2. The election 




was 



j | was not 

made before the expiration of 19 months from the priority date or, where Rule 32 applies, within the time limit under 
Rule 32.2(b)- 



The International Bureau of W1PO 
34, chemin des Colombettes 
1211 Geneva 20, Switzerland 


Authorized officer 

S. Mafia 




Facsimile No.: (41-22) 740.14.35 


Telephone No.: (41-22) 338.83.38 
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REQUEST 



The undersigned requests that the present 

international application be processed 
according to the Patent Cooperation Treaty. 



ceiving Office use only 



International Application No. 



International Filing Date 



Name of receiving Office and "PCT International Application" 



Applicant's or agent's file reference 
(if desired) (12 characters maximum) 



114481.5 RS 



Box No. I TITLE OF INVENTION 

DETERMINATION OF TEMPERATURE AND/OR tMlSSlVI" 


TY FUNCTION OF OBJECTS BY REMOTE SENSING 


Box No. II APPLICANT 


Name and address: (Family name followed by given name; for a legal entity, full official 
The address must include postal code and name of country. The country of the address indicated in this 
Box is the applicant's State (that is, country) of residence if no State of residence is indicated below.) 
ISRAEL AIRCRAFT INDUSTRIES LTD., 
TAMAM DIVISION 
P.O.Box 75 
Yahud 56100 
Israel 


| | This person is also inventor. 


Telephone No. 


Facsimile No. 


Teleprinter No. 


State (thai is, country) of nationality: 
IL 


State (that is, country) of residence: 
IL 


This person is applicant | I all designated K7\ all designated States except I I the United States I I the States indicated in 
for the purposes of: I I States l^NI the United States of America I I of America only I I the Supplemental Box 


Box No. Ill FURTHER APPLICANT(S) AND/OR (FURTHER) INVENTOR(S) 


Name and address: (Family name followed by given name; for a legal entity, full official designation. 
The address must include postal code and name of country. The country of the address indicated in this 
Box is the applicant's State (that is, country) of residence if no State of residence is indicated below.) 
KOLTUNOV, Yoseph 
10/20 Avivim Street 
Holon 58267 
Israel 


This person is: 

| | applicant only 

applicant and inventor 

| | inventor only (If this check-box 
is marked, do not fill in below.) 


State (that is, country) of nationality: 
IL 


State (that is, country) of residence: 
IL 


This person is applicant I | all designated | | all designated States except K7\ the United States I I the States indicated in 
for the purposes of: I — I States I I the United States of America IcSI of America only I I the Supplemental Box 


Further applicants and/or (further) inventors are indicated on a continuation sheet. 


Box No. IV AGENT OR COMMON REPRESENTATIVE; OR ADDRESS FOR CORRESPONDENCE 


The person identified below is hereby/has been appointed to act on behalf |^| t | — | 

oftheapplicant(s) before the competent International Authorities as: l&l a S ent I I common representative 


Name and address: (Family name followed by given name; for a legal entity, full official 
designation. The address must include postal code and name of country.) 

REINHOLD COHN AND PARTNERS 
P.O.Box 4060 
Tel-Aviv 61040 
Israel 


Telephone No. 
+972 3 7109312/3 


Facsimile No. 
+972 3 5606405 


Teleprinter No. 


I I Address for correspondence: Mark this check-box where no agent or common representative is/has been appointed and the 
I — I space above is used instead to indicate a special address to which correspondence should be sent. 
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Sheet No. 



Continuation of Box No. Ill FURTHER APPLICANTS AND/OR (FURTHER) INVENTOR(S) 



If none of the following sub-boxes is used, this sheet is not to be included in the request. 


Name and address: (Family name followed by given name; for a legal entity 1 , full official designation. 

The address must include postal code and name of country. The country of the address indicated in this 

Box is the applicant's State (that is, country) of residence if no State of residence is indicated below.) 

MAXIMOV, Alexander 

17/20 Haarazim Street 

Jerusalem 96181 

Israel 


This person is: 

| | applicant only 

fc^j applicant and inventor 

| | inventor onlv (If this check-box 
is marked, do not fill in below.) 


State (that is. country) of nationality: 
IL 


State (that is, country) of residence: 
IL 


This person is applicant I I all designated I I all designated States except K7| the United States | I the States indicated in 
for the purposes of: I I States I I the United States of America le^l of America only I I the Supplemental Box 


Name and address: (Family name followed by given name; for a legal entity, full official designation. 
The address must include postal code and name of country. The country of the address indicated in this 
Box is the applicant's State (that is, country) of residence if no State of residence is indicated below.) 

MEITIN, Igor 

7/30 Ben-Gurion Street 

Giv'at Shmuel 54018 

Israel 


This person is: 

| | applicant only 

IX! applicant and inventor 

| | inventor onlv (If this check-box 
is marked, do not fill in below.) 


State (that is, country) of nationality: 
IL 


State (that is, country) of residence: 
IL 


This person is applicant I I all designated I I all designated States except K7I the United States I I the States indicated in 
for the purposes of: I — I States I I the United States of America l^j of America only I I the Supplemental Box 


Name and address: (Family name followed by given name; for a legal entity, full official designation. 
The address must include postal code and name of country. The country of the address indicated in this 
Box is the applicant's State (that is, country) of residence if no State of residence is indicated below.) 

ALLON, Motti 
44 Rabinovitz Street 
Holon 58672 
Israel 


This person is: 

| | applicant only 

IX! applicant and inventor 

| | inventor onlv (If this check-box 
is marked, do not fill in below.) 


State (that is, country) of nationality: 
IL 


State (that is, country) of residence: 
IL 


This person is applicant I I all designated I I all designated States except K7j the United States I I the States indicated in 
for the purposes of: I I States I I the United States of America IA1 of America only I I the Supplemental Box 


Name and address: (Family name followed by given name; for a legal entity, full official designation. 
The address must include postal code and name of country. The country of the address indicated in this 
Box is the applicant's State (that is, country) of residence if no State of residence is indicated below.) 

GUTTMAN, Glen D. 

39 Sderot Chen 

Tel-Aviv 64166 

Israel 


This person is: 

| | applicant only 

applicant and inventor 

| | inventor onlv (If this check-box 
is marked, do not fill in below.) 


State (that is, country) of nationality: 
IL 


State (that is, country) of residence: 
IL 


This person is applicant I I ail designated I I all designated States except K7| the United States I I the States indicated in 
for the purposes of; I — I States I I the United States of America l^l of America only I I the Supplemental Box 


[XI Further applicants and/or (further) inventors are indicated on another continuation sheet. 
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Continuation of Box No. Hf^TIRTHER APPLICANTS AND/OR (FUR 



) INVENTOR(S) 



If none of the following sub-boxes is used, this sheet is not to be included in the request. 


Name and address: (Family name followed by given name; for a legal entity, full official designation. 
The address must include postal code and name of country. The country of the address indicated in this 
Box is the applicant's State (that is, country) of residence if no State of residence is indicated below.) 

KERSHENBAUM, Arik 
3/1 Moshe Dayan Street 
Ra'anana 43580 
Israel 


This person is: 

| | applicant only 

applicant and inventor 

| | inventor onlv (If this check-box 
is marked, do not fill in below.) 


State (that is, country) of nationality: 
IL 


State (that is, country) of residence: 
IL 


This person is applicant I I all designated I | all designated States except K7| the United States I I the States indicated in 
for the purposes of: I I States I I the United States of America k^A of America only I I the Supplemental Box 


Name and address: (Family name followed by given name; for a legal entity, full official designation. 
The address must include postal code and name of country. The country of the address indicated in this 
Box is the applicant's State (that is, country) of residence if no State of residence is indicated below.) 


This person is: 
i — i 

| j applicant only 

| | applicant and inventor 

| j inventor onlv (If this check-box 
is marked, do not fill in below.) 


State (that is, country) of nationality: 


State (that is. country) of residence: 


This person is applicant I I all designated I I all designated States except | I the United States | I the States indicated in 
for the purposes of: 1 1 States 1 1 the United States of America 1 1 of America only 1 I the Supplemental Box 


Name and address: (Family name followed by given name; for a legal entity, full official designation. 
The address must include postal code and name of country. The country of the address indicated in this 
Box is the applicant's State (that is, country) of residence if no State of residence is indicated below.) 


This person is: 

| | applicant only 

| | applicant and inventor 

j | inventor onlv (If this check-box 
is marked, do not fill in below.) 


State (that is, country) of nationality: 


State (that is, country) of residence: 


This person is applicant 1 1 all designated 1 1 all designated States except I 1 the United States | I the States indicated in 
for the purposes of: I — I States 1 1 the United States of America 1 1 of America only 1 1 the Supplemental Box 


Name and address: (Family name followed by given name; for a legal entity, full official designation. 
The address must include postal code and name of country. The country of the address indicated in this 
Box is the applicant's State (that is, country) of residence if no State of residence is indicated below.) 


This person is: 

| | applicant only 

| | applicant and inventor 

| | inventor onlv (If this check-box 
is marked, do not fill in below.) 


State (that is, country) of nationality: 


State (that is, country) of residence: 


This person is applicant I j all designated 1 1 all designated States except | 1 the United States 1 1 the States indicated in 
for the purposes of: 1 1 States 1 1 the United States of America 1 1 of America only 1 I the Supplemental Box 


| | Further applicants and/or (further) inventors are indicated on another continuation sheet. 
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Box No.V DESIGNATION 



Sheet No. . . .4. 



'ATES 



The following designations are hereby made under Rule 4.9(a) (mark the applicable check-boxes; at least one must be marked): 
Regional Patent 

AP ARIPO Patent: GH Ghana. GM Gambia, KE Kenya. LS Lesotho. MW Malawi, SD Sudan, SZ Swaziland, 
UG Uganda, ZW Zimbabwe, and any other State which is a Contracting State of the Harare Protocol and of the PCT 



IS 

181 



EA Eurasian Patent: AM Armenia, AZ Azerbaijan, BY Belarus. KG Kyrgyzstan, KZ Kazakhstan, MD Republic of 
Moldova, RU Russian Federation, TJ Tajikistan, TM Turkmenistan, and any other State which is a Contracting State of 
the Eurasian Patent Convention and of the PCT 

EP European Patent: AT Austria, BE Belgium, CH and LI Switzerland and Liechtenstein, CY Cyprus, DE Germany, DK 
Denmark, ES Spain, FI Finland, FR France, GB United Kingdom, GR Greece, IE Ireland, IT Italy, LU 
Luxembourg, MC Monaco, NL Netherlands, PT Portugal, SE Sweden, and any other State which is a Contracting State 
of the European Patent Convention and of the PCT 

Ig OA OAPI Patent: BF Burkina Faso, BJ Benin, CF Central African Republic, CG Congo, CI Cote d'lvoire, CM Cameroon, 
GA Gabon, GN Guinea, ML Mali, MR Mauritania. NE Niger, SN Senegal, TD Chad, TG Togo, and any other State 
which is a member State of OAPI and a Contracting State of the PCT (if other kind of protection or treatment desired, 
specify on dotted line) 

National Patent (if other kind of protection or treatment desired, specify on dotted line): 



81 
81 
Si 
81 
81 
81 
81 

b 

81 
81 
81 
8! 
81 
81 
81 
81 
81 
81 
81 
81 
81 
81 
81 
81 
□ 
81 
81 
81 
8! 
81 
81 
81 
81 
81 

81 
8! 
81 
81 
81 



AL Albania 

AM Armenia 

AT Austria 

AU Australia 

AZ Azerbaijan 

BA Bosnia and Herzegovina 

BB Barbados 

BG Bulgaria 

BR Brazil 

BY Belarus 

CA Canada 

CH and LI Switzerland and Liechtenstein 

CN China 

CU Cuba 

CZ Czech Republic 

DE Germany 

DK Denmark 

EE Estonia 

ES Spain 

FI Finland 

GB United Kingdom 

GE Georgia 

GH Ghana 

GM Gambia 

GW Guinea-Bissau 

HR Croatia 

HU Hungary 

ID Indonesia 

IL Israel 

IS Iceland 

JP Japan 

KE Kenya 

KG Kyrgyzstan 

KP Democratic Peop le's Repu bl ic of Korea 



LS Lesotho 

LT Lithuania 
LU Luxembourg 
LV Latvia 

MD Republic of Moldova 

MG Madagascar 

MK The former Yugoslav Republic of Macedonia 



KR 
KZ 
LC 
LK 
LR 



Republic of Korea 

Kazakhstan 

Saint Lucia 
Sri Lanka 
Liberia 



8] 
8J 
81 
81 
8! 
81 
81 

81 
81 
81 
8! 
8! 
81 
81 
8! 
8! 
81 
81 
81 
8! 
81 
81 
81 
81 
81 
81 
81 
81 
81 

81 
81 
8! 
8i 

Check-boxes reserved for designating States (for the purposes of 
a national patent) which have become party to the PCT after 
issuance of this sheet: 



MN Mongolia 

MW Malawi 

MX Mexico 

NO Norway 

NZ New Zealand 

PL Poland 

PT Portugal 

RO Romania 

RU Russian Federation 

SD Sudan 

SE Sweden 

SG Singapore 

SI Slovenia 

SK Slovakia 

SL Sierra Leone 

TJ Tajikistan 

TM Turkmenistan 

TR Turkey 

TT Trinidad and Tobago . . 

UA Ukraine , 

UG Uganda 

US United States of America 



UZ 
VN 
YU 



Uzbekistan 
Viet Nam 
Yugoslavia 



ZW Zimbabwe 



IS 

□ 



GD- Grenada , 



Precautionary Designation Statement: In addition to the designations made above, the applicant also makes under Rule 4.9(b) all 
other designations which would be permitted under the PCT except any designation(s) indicated in the Supplemental Box as being 
excluded from the scope of this statement. The applicant declares that those additional designations are subject to confirmation and 
that any designation which is not confirmed before the expiration of 15 months from the priority date is to be regarded as withdrawn 
by the applicant at the expiration of that time limit. (Confirmation of a designation consists of the filing of a notice specifying that designation 
and the payment of the designation and confirmation fees. Confirmation must reach the receiving Office within the 1 5-month time limit.) 
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Sheet No. 



Box No. VI PRIORITY C 



Filing date 
of earlier application 
(day/m onth/year) 



m 



Number 
'earlier application 



[further priority claims are indicated in the Supplemental Box. 
irlier application is: 



national application: 
country 



re^rona 



fhal application:* 
regional Office 



international application: 
receiving Office 



item(l) 

20 November 1997 

fro. VI 



122258 



IL 



item (2) 



item (3) 



The receiving Office is requested to prepare and transmit to the International Bureau a certified copy 
of the earlier application(s) (only if the earlier application was filed with the Office which for the 
purposes of the present international application is the receiving Office) identified above as item(s): 



* Where the earlier application it an AIUPO application, it is mandatory to indicate in the Supplemental Box at least one country party to the Paris Convention for the 
Protection of Industrial Property for which that earlier application was filed (Ride 4. i0(b) See Supplemental Box. 



Box No. VH INTERNATIONAL SEARCHING AUTHORITY 



Choice of International Searching Authority (ISA) 

(if two or more International Searching Authorities are 
competent to carry out the international search, indicate the 
Authority chosen; the two-letter code may be used): 

ISA/ep 



Request to use results of earlier search; reference to that search (if an earlier 
search has been carried out by or requested from the International Searching Authority): 

Date (day/month/year) Number Country (or regional Office) 



Box No. VHI CHECK LIST: LANGUAGE OF FILING 



This international application contains 
the following number of sheets : 

request 

description (excluding) 
sequence listing part) 

claims 

abstract 

drawings 

sequence listing part 
of description 

Total number of sheets 



61 



This international application is accompanied by the item(s) marked below: 





i El 


5 


2. EI 




3.D 


36 


4.D 


11 


s.n 


1 


«.□ 


8 


7.D 




«.□ 




9.D 



Figure of the drawings which y 
should accompany the abstract: T 


Language of filing of the 

international application: ENGLISH 


Box No. IX SIGNATURE OF APPLICANT OR AGENT 



Next to each signature, indicate the name of the person signing and the capacity in which the person signs (if such capacity is not 
obvious from reading the request). - 




LEVIN Elena 
Agent for REINHO 



COHN AND PARTNERS 



1. 


Date of actual receipt of the purported 
international application: 


2. Drawings: 


3. 


Corrected date of actual receipt due to later but 
timely received papers or drawings completing the 
purported international application: 




| | received: 


4. 


Date of timely receipt of the required 
corrections under PCT Article 1 1(2): 




| 1 not received: 


5. International Searching Authority tq A / 
(if two or more are competent): lo A/ 


6. | 1 Transmittal of search copy delayed 
L— ! until search fee is paid 





For International Bureau use only 



Date of receipt of the record copy 
by the International Bureau: 
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jVTENT COOPERATION TREAT 

PCT 



INTERNATIONAL SEARCH REPORT 

(PCT Article 18 and Rules 43 and 44) 



Applicant's or agent's file reference 

114481.5 RS 


FOR FURTHER see Notification of Transmittal of International Search Report 
ACTION * F ° rm PCT/ISA/220 > as we " as ' where applicable, item 5 beiow. 


International application No. 

PCT/IL 98/00568 


International filing date (day/month/year) 

20/1 1/1998 


(Earliest) Priority Date (day/month/year) 

20/11/1997 


Applicant 

ISRAEL AIRCRAFT INDUSTRIES LTD. et al . 



This International Search Report has been prepared by this International Searching Authority and is transmitted to the applicant 
according to Article 18. A copy is being transmitted to the International Bureau. 



. sheets. 



This International Search Report consists of a total of 3 

PH '* is also accompanied by a copy of each prior art document cited in this report. 



1 . Basis of the report 

a. With regard to the language, the international search was carried out on the basis of the international application in the 
language rn which it was filed, unless otherwise indicated under this item. 

the international search was carried out on the basis of a translation of the international application furnished to this 
Authority (Rule 23.1 (b)). 

b. With regard to any nucleotide and/or amino acid sequence disclosed in the international application, the international search 
was carried out on the basis of the sequence listing : 

I | contained in the international application in written form. 

I I filed together with the international application in computer readable form. 

I | furnished subsequently to this Authority in written form. 

I I furnished subsequently to this Authority in computer readble form. 

□ 

the statement that the subsequently furnished written sequence listing does not go beyond the disclosure in the 
international application as filed has been furnished 

the statement that the information recorded in computer readable form is identical to the written sequence listing has been 
furnished 

2. Certain claims were found unsearchable (See Box I). 

3 £~J Unity of invention is lacking (see Box II). 

4. With regard to the title, 

|X] the text is approved as submitted by the applicant. 

I | the text has been established by this Authority to read as follows: 



5. With regard to the abstract, 

PH the text is approved as submitted by the applicant. 

[~H the text has been established, according to Rule 38.2(b), by this Authority as it appears in Box III. The applicant may 
within one month from the date of mailing of this international search report, submit comments to this Authority. 

6. The figure of the drawings to be published with the abstract is Figure No. J 

E as suggested by the applicant. rj None of the figures. 

I | because the applicant failed to suggest a figure. 

I I because this figure better characterizes the invention. 
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INTERNATIONAL SEARCH REPORT 



A. CLASSIFICATION OF SUBJECT MATi cri 

IPC 6 G01J5/60 



International Application No 

^^IL 98/00568 



According to Intsmational Patent Classification (IPC) or to both national classification and IPC 



B. FIELDS SEARCHED 



Minimum documentation searched (classification system followed by classification symbols) 

IPC 6 G01J 



Documentation searched other than minimum documentation to the extent that such documents are included in the fields searched 



Electronic data base consulted during the international search (name of data base and, where practical, search terms used) 



C. DOCUMENTS CONSIDERED TO BE RELEVANT 



Category ° 


Citation of document, with indication, where appropriate, of the relevant passages 


Relevant to claim No. 


X 


EP 0 143 282 A (DEUTSCHE FORSCH LUFT 


1-6,21, 




RAUMFAHRT) 5 June 1985 


23, 






29-32, 






45,46 




see page 5, line 6 - page 19, line 14 






see page 22, line 19 - line 20 




X 


EP 0 420 108 A (EURATOM) 3 April 1991 


1-6, 






8-10,29, 






30 




see the whole document 




X 


US 5 377 126 A (FLIK MARKUS I ET AL) 


1-6, 




27 December 1994 


8-11,29, 






30 




see the whole document 






-/-- 





Further documents are listed in the continuation of box C. 



Patent family members are listed in annex. 



° Special categories of cited documents : 

"A" document defining the general state of the art which is not 

considered to be of particular relevance 
"E" earlier document but published on or after the international 

filing date 

"L" document which may throw doubts on priority claim(3) or 
which ts cited to establish the publication date of another 
citation or other special reason (as specified) 

"O" document referring to an oral disclosure, use, exhibition or 
other means 

"P" document published prior to the international filing date but 
later than the priority date claimed 



T" later document published after the international filing date 
or priority date and not in conflict with the application but 
cited to understand the principle or theory underlying the 
invention 

"X" document of particular relevance; the claimed invention 
cannot be considered novel or cannot be considered to 
involve an inventive step when the document is taken alone 

"Y" document of particular relevance; the claimed invention 

cannot be considered to invotve an inventive step when the 
document is combined with one or more other such docu- 
ments, such combination being obvious to a person skilled 
in the art. 

document member of the same patent family 



Date of the actual completion of the international search 

2 March 1999 


Date of mailing of the international search report 

12/03/1999 


Name and mailing address of the ISA 

European Patent Office, P.S. 5616 Patentlaan 2 
NL - 2280 HV Rijswijk 
Tel. (+31-70) 340-2040, Tx. 31 651 epo nl. 
Fax: (+31-70) 340-3016 


Authorized officer 

Ramboer, P 
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INTERNATIONAL SEARCH REPORT 



^Continuation) DOCUMENTS CONSIDERED TO BE RELEVANT 



International Application No 

^^/IL 98/00568 



Categor/ ° Citation of document, with indication, where appropriate, of the relevant passages 



Relevant to claim No. 



GB 2 131 649 A (BRITISH AEROSPACE) 

20 June 1984 

see the whole document 

WO 97 11340 A (3T TRUE TEMPERATURE 
TECHN0L0GI ;GLAZMAN EVGENY (IL)) 
27 March 1997 

see page 8, line 3 - page 9 5 line 4 
see page 9, line 15 - line 38 

KEALY P S ET AL: "SEPARATING TEMPERATURE 
AND EMISSIVITY IN THERMAL INFRARED 
MULTISPECTRAL SCANNER DATA: IMPLICATIONS 
FOR RECOVERING LAND SURFACE TEMPERATURES" 
IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE 
SENSING, 

vol. 31, no. 6, 1 November 1993, paqes 
1155-1164, XP000432082 
see the whole document 

DE 41 34 313 A (DEUTSCHE AEROSPACE) 

22 April 1993 

see the whole document 

RIICHI NAGURA ET AL: "OPTICAL SENSOR 
SYSTEM FOR JAPANESE EARTH RESOURCES 
SATELLITE-1" 

NEC RESEARCH AND DEVELOPMENT, 
vol. 33, no. 2, 1 April 1992, pages 
226-235, XP000310833 
see the whole document 

DATABASE WPI 

Section EI, Week 8839 

Derwent Publications Ltd., London, GB; 

Class S03, AN 88-277509 

XP002095191 

& SU 1 383 259 A (AEROGEOLOGIYA) 
, 23 March 1088 
see abstract 



1,12, 
50-53 



1,2,5 



1,2,12, 

29,30, 

50,51 



1,50-53 



1,4,12, 

23,29, 

45,46,53 



1,12,13, 

25,26, 

29,50-53 
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Infj^Afcjon on patent family members 



Patent document 
cited in search report 



Publication 
date 



Patent family 
member(s) 



International Application No 

^^IL 98/00568 



Publication 
date 







A 

n 




DE 


3414984 


A 


24-10-1985 


DD 


229213 


A 


30-10-1985 


JP 


1796595 


C 


28-10-1993 


OP 


5005051 


B 


21-01-1993 


JP 


60140131 


A 


25-07-1985 


US 


4974182 


A 


27-11-1990 


DD 


229212 


A 


30-10-1985 



EP 0143282 



05-06-1985 



EP 0420108 A 03-04-1991 



LU 


87595 


A 


07- 


-05- 


-1991 


AT 


107021 


T 


15- 


-06- 


-1994 


AU 


639029 


B 


15- 


-07- 


1993 


AU 


6418590 


A 


18- 


-04- 


-1991 


CA 
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(57) Abstract 



A method and a system for accurate determining temperature and/or emissivity of an object by remote sensing, comprising acquiring 
electromagnetic radiation data from the object in N spectral bands, deriving temperature and/or emissivity values from said data using a 
system of equations, each describing the electromagnetic radiation that is emitted from the object in the particular band as a function of 
temperature and emissivity of the object, wherein the emissivity depends on temperature and wavelength. Based on the above concept 
of temperature and emissivity determination, there are' also developed a method and a system for the object's pattern recognition, being 
applicable to a number of purposes. 
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Determination of temperature and/or emissivity function of objects by 

remote sensing 

5 

FIELD OF THE INVENTION 

The present invention relates to a method and a system for 
remotely determining temperature and emissivity parameters of objects by 
multispectral measurements. 

10 BACKGROUND OF THE INVENTION 

Various approaches have been used, for example, for passive 
remote sensing of the ground. Thus, maps of surface brightness, temperature 
and emissivity have been used for investigating geological surface properties 
(Kahle et aL 1980. Applied Optics, 19. 2279), whereas maps of thermal 

15 inertia have been used to infer subsurface properties of soil (Price, 1977, 
Journal of Geophysical Research, 82, 2582). 

The detection of underground structures requires penetration of the 
ground and is therefore accomplished by active sensing techniques, such as 
radar (Blake. 1993, "Ground-Penetration Radar Developed by Sweden", 

20 International Defense Review, 3, 193; von Maydell et aL 1987, U.S. Patent 
4.675,677), or combined passive sensing and radar (Clark et aL SPIE - The 
International Society for Optical Engineering; 1942, 178). 

It is acknowledged in the prior art. that emissivity of various 
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materials may noticeably depend on wavelength of the emitted radiation and 
temperature of the sample. 

For example. US 4,659,234 describes a method of emissivity error 
correcting for radiation thermometer, where temperature is determined based 
5 on two measurements of radiated infrared energy at two closely adjoining 
wavelengths, and the result is corrected using one single measurement. While 
the method is aimed to solve a problem of determining the randomly changing 
value of emissivity of heated metal objects, it does not take into consideration 

nnv ripnpnrlfnrp r\f p.mi<;<sivitv nn thp rnriiatinn wnvplpncrth anH mnrh Ircc - nn 

. r - 'J ~ CT-" t 

io temperature of the object. Such an approach renders the method inaccurate. 

US 5,132.922 describes an emissivity independent 
multi-wavelength pyrometer operating according to a method using a 
least-squares-based multiwavelength pyrometry technique and a theoretical 
function for the dependence of the radiance on the wavelength. However, in 

is this method the emissivity /wavelength function is considered to be the same 
in all spectral bands, which either leads to a significant reduction in the 
accuracy of the determined temperature, especially for heated non-transparent 
materials, or to the use of narrow-band filters. The latter poses a problem for 
infra-red measurements in real time, since narrow band filters provide 

20 insufficient energy for further processing. 

Moreover, all the described above methods make an assumption 
that emissivity of an object is independent of its temperature, since none of 
the methods appear to suggest any way of considering this dependency for 
determining temperatures. Such an assumption inadvertently leads to essential 

25 errors in the temperature determination for at least non-transparent heated 
material. 

In many branches of industry, for example in metal lurgy or in 
manufacture of semiconductor devices samples and materials are exposed to 
various operations while being in the heated condition. 
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Determining temperature of such objects during the manufacturing 
process presents a very important task, which is rather difficult due to the 
following factors: 

- it is impossible to perform direct temperature measurements, so 
5 the temperature determining is usually performed via the remotely acquiring 

data on the infrared radiation of the sample and further processing this data, 

- in the mathematical processing of the remotely obtained data 
many parameters intrinsic to the material of interest must be used (such as its 
thermal emissivity. etc.)., which are usually unknown and cannot be directly 

io measured. 

Another important task is detection (usually in the dark) and 
recognition of various objects having different temperatures and/or 
emissivities. Difficulties of such a process are essentially the same as outlined 
above. 

15 

SUMMARY OF THE INVENTION 

It is therefore the purpose of the present invention to provide an 
accurate method of remotely determining temperature and/or emissivity 
function of objects. It is to be noted, that the problem of detemiining the 
20 emissivity function of real objects is one of important problems in many 
industrial applications. 

Another objective of the invention is designing a system for 
implementing the above-mentioned method. 

Yet another objective of the present invention is to provide a 
25 method and a system being capable of detecting and. where applicable, 
identifying objects on the basis of remotely determined temperature and 
emissivity distributions. 

A specific objective of the present invention will be to provide a 
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real time recognition of objects (usually in the dark) from a movable platform, 
using the inventive method. 

The term "object" is used herein in a broad sense to denote any 
kind of structure, crop or articles in a solid, liquid or gaseous state on the 
5 ground, natural and artificial underground structure, as well as objects located 
above the ground (such as solid samples, oil spots on water, clouds, overhead 
electric wires, etc.). 

For the sake of simplicity the term "object detection" will be used 
in the following description and claims to denote the detection of an object or 
10 objects or of some properties thereof. 

The above main purpose can be achieved by providing a method 
of remotely determining temperature of an object, comprising: 

- acquiring electromagnetic radiation data from the object in N spectral 
bands, 

15 - deriving temperature value from said data using a system of equations, 

each describing the electromagnetic radiation that is emitted from the object 
in said band as a function of temperature and emissivity of the object: 
the method being characterized in that, 

- said electromagnetic radiation is measured M times, each time in N 
20 said spectral bands, thus obtaining N*M readings, 

- a system of N*M equations is formed, each describing the 
electromagnetic radiation that is emitted from the object in said band as a 
function of temperature T and emissivity e of the object, wherein the 
emissivity £ is described by a model incorporating at least one function, of at 

25 least two arguments being temperature T and wavelength a.; said system 
comprising N*M readings of the electromagnetic radiation and at least M+P 
unknowns which include M values of the object's temperature corresponding 
to said respective M measurements of the electromagnetic radiation and P 
parameters, wherein P > 2; said parameters belonging to said model of 
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emissivity; 

solving said system of N*M equations, thus deriving temperature 
of the object at any one of the M moments of said electromagnetic radiation 
measurement, wherein M is a natural number of the set [L 2, 3. . .]. 
5 According to one particular version of the above method, the 

mentioned model is selected to comprise one function s = f (A*T) taking 
constant values within respective spectral bands, and therefore P is considered 
equal to N, 

therebv said svstem of N*M equations is formed to comDrise N*M readings 

- * A 

io of the electromagnetic radiation and N+M unknowns which include M values 
of the object's temperature corresponding to said respective M measurements 
of the electromagnetic radiation, and N values of emissivity for said 
respective N spectral bands. 

The above-defined particular method is based on the following 

15 assumptions: 

1 . Emissivity of the object in each one of the N spectral bands is a 
predetermined function of temperature having parameters which are constant 
during all the M measurements of the electromagnetic radiation performed in 
this band. 

20 2. During each specific measurement of the electromagnetic 

radiation, temperature of the sample is assumed to be constant (i.e. its change 
is considered negligible), regardless the spectral band where such specific 
measurement is performed. 

The assumption 2 is best fulfilled if each measurement of the 

25 electromagnetic radiation is performed in all N spectral bands substantially 
simultaneously. It is especially important for accurate measurements 
concerning objects which undergo either heating or cooling. For an object or 
a plurality of objects which have rather constant or slowly changing 
temperatures, the measurements in N spectral bands may be accomplished in 
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sequence, at a maximally available speed. 

In practice, i.e. when digital computing technology is used for 
controlling the process of measurements and processing the obtained data- 
each of said M multispectral measurements of the electromagnetic radiation 
5 (i.e. each of said M*N readings thereof) is digitized and stored in a computer 
memory for further computerized processing, said processing comprising 
substituting the digitized M*N readings of the electromagnetic radiation into 
said M*N equations, and solving thereof. 

Preferably, the above measurements are to be provided in the 
10 infrared spectrum. 

According to one specific version of the method said equations are 
integral ones, for example Fredholm equations relating to emissivity. 

It should be emphasized, that the above-defined method easily 
provides sufficiently accurate temperature values both in cases when 
15 emissivity of the object essentially varies from one spectral band to another 
(either due to dispersion between these bands, or due to specific physical 
properties of the object), and in cases when the emissivity is rather uniform 
throughout the spectrum. 

In addition, a succession of temperature values obtained by the 
20 method can be used for real-time monitoring and controlling of the object's 
thermal condition, for example during a technological process to which the 
object is subjected. 

It is well known, that at high temperatures (i.e. those between the 
room temperature and 1,500°C) emissivity of some real materials, such as 
25 semiconductor samples, varies not only from one spectral band to another, but 
also depends on the temperature of the sample. One of the reasons for such an 
effect is degrading of the sample's transparency with the temperature increase, 
that consequently affects emissivity of the sample. In other words, the 
simplest case of assumption 1 mentioned above is not applicable to such 
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objects, and. when temperature of heated non-transparent objects must be 
determined, the method where only N+M unknowns are sought for will not be 
sufficiently accurate. It means, that if the sample emissivity is considered to 
be independent from temperature while the sample undergoes any thermal 
5 treatment, this assumption will be a source of essential errors in calculations 
of the temperature. 

Having assumed, that emissivity of the heated sample is a known 
mathematical function of temperature and wavelength with unknown 
parameters, the inventors suggested to introduce such parameters as 
io unknowns into the above-mentioned system of equations, in order to obtain 
values of the parameters by solving the system. 

More particularly, each of said M*N equations describes the 
electromagnetic radiation that is emitted from the object in one of said bands 
as a runction of the temperature and the respective emissivity of the object, 
!5 the emissivity, in turn, being a known function of temperature, wavelength 
and a known number of unknown parameters; said parameters being 
introduced as additional unknowns into said system of equations for further 
obtaining their values upon solving the system. 

In one example, the mentioned known function may be 
20 represented as a sum of a known number of terms of two other functions, a 
first of which being dependent only on temperature and a known number of 
unknown parameters, and a second - only on wavelength and another known 
number of other unknown parameters. 

According to another example, this known runction may be 
25 represented by a polynomial expansion having K unknown coefficients, said 
coefficients being introduced into the system of equations as unknown 
parameters to be determined among said unknowns upon solving said system. 

If the emissivity variable in each of the above mentioned equations 
is represented as a polynomial function having K unknown coefficients, the 
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s\ stem of the equations will overcome some transformations. Namely, though 
the number of equations N*M is unchanged, the number of unknowns in the 
system will increase and become equal to M + K*R where K is the number 
of polynomial coefficients of the function emissivity/temperature in each of 
5 the N spectral bands. However, if the numbers of the polynomial coefficients 
are different in different spectral bands (which case is more general and more 
correct from the physical point of view), the number of unknowns in the 
system will be M + (K Ai + K A2 --- + K^n), 

where K A) . ..K a n are numbers of the polynomial coefficients of the 

io emissivity/temperature Function in each of the spectral bands, respectively. 

In practical applications, where the electromagnetic measurements 
are performed and processed by the aid of computer, the number of 
electromagnetic radiation measurements M and the number of spectral bands 
N may be as great as is wished. For example, in industrial conditions during a 

15 relevant technological process up to 2,000 measurements may be performed 
in a short period of time, each measurement taking from 1 to 10 msec. On the 
other hand, the number K of polynomials in the function emissivity/temperat- 
ure usually does not exceed 4. Owing to that the system of such equations is 
solvable, since 

20 N*M»M + K*R or 

N*M » M + (K A , + K A2 - + K AN ). 
The most effective application of the above described modified 
method is for non-contact temperature measurement of heated non-transpar- 
ent objects, for example of semiconductor materials and samples during their 

25 manufacture. 

It should be emphasized, that, since values (parameters) of the 
emissivity function are solutions of said system of M*N equations, as well as 
the temperature values are, the disclosed above method of remotely 
measuring temperature of an object constitutes also a method of determining 
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the emissivity function of the object. 

In other words, there is provided a method of remotely 
determining emissivity function of an object, comprising 

- acquiring electromagnetic radiation data from the object in N spectral 
5 bands in the infrared spectrum by measuring it M times, each time in N said 

spectral bands, thus obtaining N*M readings, 

- forming a system of N*M equations, each describing the electromag- 
netic radiation that is emitted from the object in one of said bands as a 
function of temperature T and emissivifv c of the ohi 

to emissivity e is described by a model incorporating at least one function, of at 
least two arguments being temperature T and wavelength a.; said system 
comprising N*M readings of the electromagnetic radiation and at least M+P 
unknowns including M values of the object's temperature corresponding to 
said respective M measurements of the electromagnetic radiation, and P 

15 parameters, wherein P > 2; said parameters belonging to said model of the 
emissivity: 

- solving said system of N*M equations to derive said P parameters, 
thereby reconstructing the emissivity function of the object. 

All the above-described modifications of the method for 
20 determining temperature are applicable to the method for determining 
emissivity, mutatis mutandis. 

It should be appreciated, that the above described methods for 
remotely determining temperature and/or emissivity function of an object may 
be utilized for deriving temperature and/or emissivity descriptive maps of the 
25 object. 

To the above-mentioned purpose, the object is to be divided into a 
plurality of representative segments, and each segment is to be investigated 
according to the above-suggested technique. The processed data on the 
temperature and/or emissivity function associated with each specific segment 
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of the object can then be used for forming a corresponding pixel of the 
respective descriptive map. 

As has been mentioned earlier, ''the object'' is defined in a broad 
sense, so it can be understood as a plurality of so-called sub-objects having 
5 different temperatures and/or emissivities and dislocated in a selected region 
which limits the object in space (for example, the sub-objects being such as 
the ground, the space above including one or more articles, the surrounding 
atmosphere in the selected region). Consequently, when the object is divided 
into a plurality of representative segments, in most cases the sub-objects will 
io not coincide with the representative segments, for example one sub-object 
may be covered by a number of segments and should be thereby represented 
by a number of corresponding pixels. 

Further, the above-described method can serve for the object 
detection on a background of the ground and/or the surrounding atmosphere 
15 in a selected region. For this purpose, the method may additionally comprise 
the following steps: 

- acquiring meteorological data indicative of climatic conditions of the 
selected region on the ground and/or the surrounding atmosphere: 

- recording the acquired electromagnetic radiation and meteorological 
20 data by means of suitable data recording means; 

- storing the so-recorded data in suitable data storage means and 
deriving from the so-stored data, descriptive maps of the selected 
region: and 

- feeding data of said descriptive maps into pattern recognition 
25 processor means to classify the pixels of said maps of said selected 

region, whereby the desired object detection is achieved. 
By one embodiment of the invention said meteoroloaical data are 
obtained from meteorological stations. 

Typically the said descriptive maps are temperature, emissivity and 
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albedo maps. Albedo is a parameter which can be defined as follows: 

where. r v is the radiation measured by the radiometer in the visible and near 
5 infra red regions; Tv , is the atmospheric transmissivity in the visible and near 
infra red regions; To , is the optical transmissivity of the radiometer including 
the transmissivity of the spectral filter in the visible and near infra red regions: 
and W s is the radiation from the sky. 

According to one embodiment of the invention applicable in cases 
10 where the objects detection is based on pronounced temperature and/or 
emissivity differential between the representative segments of the object, the 
object detection according to the invention is based on a single scanning pass 
of the selected region (i.e. M=l). This embodiment of the invention may be 
applicable, for example, to the detection of overhead high voltage wires, 
15 defects in the wires, connectors and insulators, or diseases in crops. In the 
performance of this embodiment, the descriptive maps describe temperature 
and emissivity and the pattern recognition processor means classifies the 
pixels according to temperature and/or emissivity difference. Typically, the 
emissivity map is used by the partem recognition means in the classification 
20 process. In the case of overhead high voltage systems for example, any such 
local temperature increase denotes a possible defect in the inquiry sector 
thereof. 

If desired, the emissivity map can be used to remove background 
noise, hence enhancing the temperature differences of the object under 
25 investigation. This is done by using the emissivity map to identify those 
pixels belonging to the object under investigation. Any temperature 
differences belonging to the background are then eliminated. 
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According to another embodiment of the invention the object 
detection is based on multiple scanning passes at desired times (M=2. 3. 4...) 
and said descriptive maps describe temperature, emissivity and albedo, 
\\ herefrom a set of characteristic maps of the selected region are derived by 
5 incorporating suitable reference data. 

Typically an optimal set of characteristic maps of the selected 
region is chosen and data thereof are fed into the pattern recognition 
processor means in order to classify the pixels of said optimal set of 
characteristic maps of the selected region, whereby the desired object 
io detection is achieved. The set of characteristic maps may comprise maps of 
thermal inertia, thermal flux, coefficient of heat transfer and coefficient of 
mass transfer. 

The scanning of the selected region for measuring electromagnetic 
radiation in a number of spectral bands may be effected either by remote 

15 scanning or by close range scanning. 

In general, the pattern recognition processor means first defines 
the number of classes that the selected region is to be divided into. It then 
divides the selected region into the defined number of classes by statistically 
associating each pixel of the characteristic maps with a given class. Although 

20 the classes are not necessarily identified with areas having different physical 
properties- e.g.. areas of water, rocks, vegetation, sand, soil, etc., there should 
be some correlation between the two for the process of associating the pixels 
to each class to be successful. For a given selected region there exist 
preliminarily determined universal reference data that are required in order to 

25 evaluate the characteristic maps. The values of the reference data are obtained 
in controlled experiments as described below. It has been found that sufficient 
accuracy can be obtained in the detection process by using an average set of 
universal constants. This set is derived by taking the average value of the set 
obtained for ground relative humidities less that 12% and that obtained for 
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ground relative humidities greater or equal to 12%. The term "universal 
reference data" indicates that, to a good approximation, the reference data are 
suitable for all kinds of backgrounds and objects. In the advent that greater 
accuracy, with respect to the ground humidity, is desired in the detection 
5 process, the universal sets from which the average set was derived can be 
used, depending on the humidity of the ground. 

When the object detection according to the invention aims at 
detecting below ground level characteristics of the selected region, e.g. at 
finding underground artificial structures, geological structures holding 
io deposited minerals, etc.. the underground structure of interest will then show 
up as an anomalous pattern in the pattern recognition process, due to the fact 
that the heating and cooling processes of the ground are different in the region 
of the underground structure as compared to the surrounding region which is 
absent of such underground structures. In other words, the pattern recognition 

15 means will not be able to associate those pixels identified with the 
underground structure with one of the defined classes. 

In general then, object detection entails pixel classification, in 
which a number of classes is defined and the pixels in the selected region are 
statistically associated with the classes. To this end a quality function is 

20 defined and the pattern recognition process entails detennining that set of 
characteristic maps that maximizes the quality function. Anomalous patterns 
are those in which the pixels cannot be associated with a given class in the 
selected region. Where possible, it is preferable that the pixel classes 
correspond to areas having well known physical properties. 

25 Although the reference data can be considered as universal, more 

accurate values can be used if desired. Values for the reference data for 
different backgrounds or objects can be obtained by measuring the radiation 
from a reference sample of known emissivity, of the background or object of 
interest, at predetermined times and at the same time directly measuring the 
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temperature thereof. The reference data can then be derived from the 
measured data by utilizing an energy balance equation. 

By one embodiment of the invention the remote scanning of the 
selected ground region is performed with at least five sensors operating, 
5 respectively, in the 3-4 ji, 4-5 u., 8-10 u, 10-12 u. and the near infrared and 
visible spectral bands, and there is provided at least one further sensor for 
measuring electromagnetic radiation from the space above, operating in the 
0.4-2 jli spectral band. 

Quite generally, the characteristic maps of the selected region 
10 include maps of thermal inertia, thermal conduction flux, coefficient of heat 
transfer, coefficient of mass transfer, and the mean, variance and texture maps 
derived therefrom. However, if the relative humidity of the ground in the 
scanned region is greater that 30%, these maps cannot be computed. 

Preferably, the said classifying processor means in which the 
15 pixels of the characteristic maps of the selected region are classified into 
classes, makes use of a poly-Gauss expansion to describe the statistical 
distribution of the pixels and Bayes 1 method to sort the pixels into classes. 

For the said evaluation of the sets of characteristic maps for the 
purpose of selecting therefrom an optimal set. the said pattern recognition 
20 means comprises quality factor evaluation means and the pattern recognition 
process entails determining a set of characteristic maps for which the quality 
factor is the highest, thereby giving rise to the best pixel classification. 

The pattern recognition means can also include additional 
geometrical classification means for determining the shape of the detected 
25 object. 

Still further, the method can be specifically adapted for real time 
temperature/emissivity mapping of a plurality of selected regions to be 
observed, in said method the step of acquiring said electromagnetic radiation 
data is provided in real time by equipment mounted on a movable platform, 
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while displacing along said plurality of selected regions; wherein the steps of 
recording and storing the acquired electromagnetic radiation and meteorologi- 
cal data and further deriving from the so-stored data, descriptive maps of the 
selected region are registered (i.e. aligned) with the platform displacement. 
5 It is understood, that the above method may include automatic 

recognition of objects positioned along or on the course of the moving 
platform, such as a vehicle. Still further, the method may comprise steps of 
initiating an alarm and. optionally, of automatic adjusting the course of the 
\ chicle whenever an obstacle is recognized. 

10 

In accordance with a further aspect of the invention, there is 
provided a system for remotely determining temperature and/or emissivity 
function of an object, said system comprising: 

scanning means for acquiring electromagnetic radiation data from 
15 the object M times, each time in N spectral bands, thus obtaining N*M 
readings of the electromagnetic radiation, wherein M is a natural number of 
the set [1,2, 3 -..]; 

means for recording and storage said N*M readings in the digital 

form; 

20 computational means for processing said N*M readings of the 

electromagnetic radiation as a system of N*M equations comprising M 
unknown values of temperature and at least P unknown parameters, wherein 
P > 2; said parameters belonging to a model of the objects emissivity. the 
model comprising at least one function of at least two arguments being 

25 temperature T and wavelength X; 

display means for displaying one or more of said unknowns 
obtained upon solving said system of equations, or derivatives therefrom. 

The above system for remotely determining temperature and/or 
emissivity function of an object can serve for deriving temperature and/or 
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emissivity descriptive maps of the object, if said scanning means are capable 
of acquiring the electromagnetic radiation from a plurality of representative 
segments of the object, and said computational means are capable of 
processing the acquired data, for deriving and displaying the temperature 
5 and or emissivity descriptive maps of the object by the display means. 

The scanning means may comprise a plurality of sensors capable 
of simultaneously measuring electromagnetic radiation in a number of 
spectral bands. 

In another embodiment, the scanning means may comprise a 
io spatial optical filter divided into a number of sites for simultaneous acquiring 
the electromagnetic radiation from said plurality of respectively located 
segments of the object; each said site having N zones responsible for splitting 
the radiation acquired by the site into N spectral bands, thereby acquiring the 
electromagnetic radiation in all N spectral bands simultaneously. The splitting 
15 can be provided by chromatic beam splitters. 

According to an alternative embodiment of the invention, the 
scanning means comprise a rotatable optical filter for acquiring 
electromagnetic radiation in said N spectral bands in sequence, and at least 
one infrared sensor. 

20 If desired, said sensors are capable of remote scanning. 

Alternatively, the sensors may be designed for close range scanning. 

By another embodiment of the invention the scanning means 
comprises at least five ground oriented sensors operating, respectively, in the 
3-4u. 4-5 lu 8-10 u, 10-12 p. and the near infrared and visible spectral bands 
25 and at least one space oriented sensor operating in the 0.4-2 u spectral band. 

In the preferred embodiment, the system is capable of object 
detecting in a selected region having a background including the ground 
and. or the surrounding atmosphere, further comprising: 

- means for acquiring meteorological data indicative of climatic 
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conditions of the selected region and the surroundings; 

- means for recording said meteorological data; 

- data storage means for storing the recorded meteorological data 
(possibly the same means as for storing the electromagnetic radiation data); 

5 - computational means for deriving from the stored 

metheorological data and the data on the electromagnetic radiation, 
descriptive maps of the selected region; and 

- pattern recognition processor means for classifying the pixels of 

said maps. 

io By one embodiment of the invention said means for acquiring 

meteorological data are associated with meteorological stations. 

In accordance with a preferred embodiment of the invention said 
computational means are capable of deriving from said descriptive maps a set 
of characteristic maps. In such an embodiment the said processor means 

15 preferably comprise optimizing means for selecting an optimal set of 
characteristic maps which is subjected to pixel classification. The pattern 
recognition processor means can comprise pixel classification means in which 
a number of classes is defined whereby the classified pixels are statistically 
associated with the classes. For example, the pattern recognition processor 

20 means can function based on a poly-Gauss expansion for determining the 
statistical distribution of the pixels, and means based on Bayes's method for 
sorting pixels into classes. The pattern recognition processor means may also 
comprise quality factor evaluation means. 

According to one specific embodiment of the invention, the 

25 pattern recognition processor means includes additional geometrical 
classification means for determining the shape of the detected object. 

In a specific embodiment of the invention, the remote sensing 
equipment and the data recording means are positioned on any movable 
platform (for example, are airborne), whereas the data storage means. 
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computational means, classification means and pattern recognition means are 
on the ground. Alternatively, in some embodiments the remote sensing 
equipment, data recording means, data storage means, computational means, 
classification means and pattern recognition means may all be positioned on a 
5 movable platform (or, in a particular case, be airborne). 

In the above embodiments, said computational means are adapted 
to stabilize and register the recorded and stored data of the acquired 
electromagnetic radiation with the platform's movement to obtain a real time 
mapping and pattern recognition. 
io It should be appreciated, that pattern recognition processor means 

can be capable of automatic recognition of objects positioned along or on the 
course of the moving platform, and if necessary, of initiating an alarm and 
automatic adjusting the course of the vehicle whenever an obstacle is 
recognized. 

15 It goes without saying, that in the embodiments where the system 

enables the descriptive maps and the pattern recognition to be accomplished, 
the display means are capable of displaying the appropriate output data. 

BRIEF DESCRIPTION OF THE DRAWINGS: 

20 In order to understand the invention and to see how it may be carried 

out in practice, a preferred embodiment will now be described, by way of 
n on- limiting example only, with reference to the accompanying drawings, in 
which: 

Figs, la and lb present a flow-chart of the algorithm illustrating a 
25 method of non-contact temperature measurement according to the invention. 

Fig. 2 is a flow-chart of a modified method based on that illustrated in 
Fig. lb. 

Fig. 3 is a block-diagram of a system implementing the inventive 
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method. 

Figs. 4a and 4b are schematic plan views of two embodiments of the 
filter shown in Fig. 3. 

Fig. 5 is a block diagram of the process for determining the character- 
5 istic maps of the region selected for scanning. 

Fis. 6 is a block diagram of the pattern recognition process. 

Fig. 7 is a block diagram of a system according to the invention for 
performing non-contact temperature/em issivity measurements and pattern 
recognition. 

io DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

Before relating to the drawings, some basic terms and indications 
of constants and variables are to be introduced, which will be used further 
below in the examples of the inventive method. 

In the case, when the method is applied to detecting and 

15 recognizing objects on or above the Earth's surface, a multi-spectral scanning 
means (e.g. radiometer) scans the area of the surface under investigation at. at 
least one pixel at a time. The radiation received from each pixel by the 
radiometer, and by a pyranometer which receives radiation from the sky, 
serves as input for processing means which computes descriptive and 

20 characteristic maps of the selected region. The processing software is based 
on the Planck "s equation: 

IV = fT A (A)T n (A)£(A,T)B(A.,T)<U+ jr A (A)T 0 U)(\-€(A 9 T))J{^T s )clA (1) 

where the integral is over the spectral region of the sensor under consideration 
25 and where: 

W is the radiation received by the radiometer, 

r A ( A ) is the spectral atmospheric transm issivity which can be obtained from 
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standard sources (Wolfe and Zissis, Editors, 1993. The Infrared Handbook - 
revised edition. Environmental Research Institute of Michigan. U.S.A.: Selby 
et aL 1978. Atmospheric Transmittance/Radiance: Computer Code 
LOWTRAN 4, report no. AFGL-TR-78-0053. Optical Physics Division, 
5 Project 7670. Air Force Geophysics Laboratory, Hanscom AFB, MA. USA). 
t.,( /.) is the spectral optical transmissivity of the radiometer including the 
transmissivity of the spectral filter, which can be measured by standard 
methods. 

£{}^T) is the eniissiviLy of a pixel in the selected region. 
10 T is the temperature of a pixel in the selected region, 

B( a.T) is the spectral radiance of a black body at surface temperature 7\ 

(for the definition of B( IJ) see - Bramson, 1968, Infrared Radiation - A 

Handbook for Applications, Plenum Press) 

7\ is the sky temperature, and 
15 Ms..TJ is the sky spectral radiance which is given by B(?.,t,), where a 

possible choice for t 5 is T s = T a - 10 °K . and where To is the temperature of 

the atmosphere in the vicinity of the selected region. 

The method will be further explained for a simplified particular 
20 example, e.g. for an object consisting of materials homogenous from the point 
of view of radiation parameters. 

Radiation of the object includes three components: 
W=W1 + W2+ W3 
Wl = e *W B 

25 W2=R *Y 

W3 - tr *G 

where: 

W - total radiation of the object received by the radiometer 
Wl - radiation emitted by the object 
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W2 - radiation reflected by the object 
W3 - radiation transmitted by the object through its body 
£ - emissiviry coefficient of the object 
R - reflectance coefficient 
5 tr - transmittance coefficient 

(l=tr + R + z ) 

Wb - radiation of the black body having temperature equal to that 
of the object 

G - radiation existing at the back side of the object 
10 Y - background radiation arriving to the object. 

For a non- transparent object (having tr = 0\ which is heated to a 
relatively high temperature with respect to the ambient one. both W2 and W3 
are equal 0. 

In this case radiation W 1 in the spectral band Ai = [Xu X*] can be 
15 represented by a formula (2) derived from the general (Planck's) equation (1): 

rr = j £ {A.T)*BU,T)*dA (2) 

i.e. taking into consideration that T ^ and x o are equal / and all members 
concerning atmosphere and sky are suppressed. 

For determining intensity of infrared radiation emitted from the object, a 
20 reference body can be selected, for which there are known its radia- 
tion/temperature dependencies in each predetermined spectral band. 

A radiometer can detect radiation of the object in several spectral bands. 
In a specific spectral band, the received radiation V may be represented by a 
Fredholm equation (3): 

A 

I" = j<p'(A)*s(A,r)*B(^T)*dA (3) 
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where 

i\){X)- function dependent on the radiometer's parameters for specific band "i" 
and on transmitting properties of the atmosphere 

If the radiation measurements are performed M times (l...j...M), and 
5 each time in N spectral bands (1...L.N), simultaneously (see blocks 1. 2. 3 in 
Fig. la), the following system of (N*M) equations may be built and solved: 

l-'^V; = \<p\A)*e{Z.T l )*B{X,T J )*dl (4) 

A, 

Solution of the above system of equations will comprise M values 
of temperature (Ti..Tj..T M ) and, (in this simplified case), N values of 
10 emissivities (si,-. si,..£n) as can be seen from both the equations (4) and Figs 
la, lb ( blocks 4 ? 5, and 6). 

Fig. 2 illustrates another final portion of the flow-chart, which, 
together with the flow-chart diagram shown in Fig. la will form a general 
algorithm for cases when the object's emissivity depends both on wavelength 
15 and temperature. The system of equations processed in block 7 of the general 
algorithm is more complex than system (4) (see also block 5 in Fig. lb), 
namely: in each of the N*M equations dependence of the emissivity variable 
"£** from temperature of the sample is taken into consideration. 

For example, emissivity function of the heated object in spectral 
20 channel "i" can be represented as follows: 

£U.T) = f^El(A)*S' k (T) (5) 

where: Sk'(T) - is a function reflecting dependency of the emissivity on 
temperature: 

Ek' (a.)- is a function reflecting dependency of the emissivity on 
25 wavelength- 

K - is the preselected number of terms S k '(T)* Ek ( X) 
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Behavior of the functions Sk'(T) and is assumed to be known ( i.e. the 
functions are selected in advance); parameters of the function E k * ( X) are also 
assumed to be known. 

Based on the above, the following system of N*M equations may be 
5 built and solved: 



Solving the system (6) permits derivation of M values of the sample 
10 temperature, and N*K parameters for each of the functions Sk'(Tj). 

In yet another example, the above mentioned emissivity function Sk'(Tj) 
may be expanded into a polynomial series, and the equation (5) will acquire 
the following form: 



(6) 



where: 




(7) 



is integral of the preselected functions E, B and <p. 



(8) 



15 where: 



b'kq = parameters of function Sk'(Tj) in the "r spectral band; 



P'q(Ti) = polynomial of a pre-selected form and with pre-selected 
coefficients; 



(b' kq * P^Tj)) forms Sk'dj); 



20 



0 - number of polynomial coefficients. 



The system of equations in this case will be as follows: 



(9) 



where: 
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i = ! LN} ; j = {1,M}, k ={ LK}, q = { 1,Q}. 

Solution of the above system of N*M equations will give M unknown 
values of temperature, andN*K*Q parameters b ; kq of the function Sk^Tj) (for 
the case when in each spectral band there is one and the same number of the 
5 parameters). 

It should be emphasized that, knowing the above mentioned parameters 
and the functions Ek'(A.), one can exactly define the wavelength and 
temperature dependent function of emissivity. 

The above system of equations is linear for parameters biJ. and 

— i - 

10 non-linear for values of temperature Tj. 

Accuracy of the temperature calculations depends on the following 
parameters: 

number N of the spectral bands, and on specific selection of the 

bands. 

15 - number M of the "temperature" (i.e. electromagnetic radiation) 

measurements, 

quality of synchronizing of the "temperature" measurements in the 
N spectral bands, i.e. performing thereof (in the best case) substantially 
simultaneously, 

20 - technical parameters of radiometers and accuracy of the radiomet- 

ric measurements. 

accuracy of the mathematical model. 
Based on the above description and examples, the method of remotely 
measuring temperature of an object suggested in the present application 
25 represents also a method of determining emissivity function of the object. 

Fig. 3 illustrates a block-diagram of a system capable for 
implementing the inventive method. An object 10 (for example, a heated 
object) emits infrared radiation on to an optics 12 via which it is transmitted 
to a multichannel spatial filter 14. The multichannel spatial filter 14 performs 
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simultaneous splitting of spectral lines into N spectral bands. In Fig. 4a there 
is shown a schematic plan view 28 of the multichannel spatial filter which is 
designed to split the spectrum into N=12 spectral bands (channels) by 12 
zones capable to transmit radiation of different wavelengths. Practical 
5 implementation of the filter is not discussed in the frame of the present patent 
application, intensities of the obtained N radiation portions are then detected 
by a detector 1 6 adapted for multichannel detection. "N" analogous readings 
of the electromagnetic radiation obtained by the detector 16 are converted into 
the digital form and transmitted to a processing unit 1 8 where they are stored 

i o for further processing. 

All the above-described operations are synchronized by a control 
unit 20 of the processing unit 18, which generates M successive synchro- 
nizing signals. Each of these M signals simultaneously activates optics 12, 
multichannel filter 14, detector 16 and processing unit 18 for performing and 

1 5 recording one of the M measurements of the infrared radiation being effected 
simultaneously in N spectral bands. M*N obtained digitized readings of the 
infrared radiation emitted from the object 10 are processed in the unit 18 as a 
system of equations having at least M+P unknowns comprising M values of 
temperature of the object and at least P parameters of its emissivity model. 

20 Values of the temperature ( and/or of the emissivity), being solutions of these 
equations, are displayed by a display 22. If desired, the system can be adapted 
for deriving a temperature descriptive map of the object, as well an emissivity 
descriptive map thereof. To this end, optional thermal mapping unit 24 and 
an emissivity mapping unit are connected to the processing unit 18. Fig. 4b 

25 shows a filter structure suitable for making the mapping operation. 

Fig. 4b illustrates a plan view of a spatial filter 30 divided into a 
number of sites 32 for acquiring infrared radiation from a number of 
respectively located segments of the object (not shown). Each site has N=4 
zones responsible for splitting the infrared radiation acquired by the site into 
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four spectral bands. Plurality of the split radiation portions from all sites of 
the filter 30 may be then transmitted to a detector similar to that shown as 16 
in Fig. 3, and further to a similar processing unit. The processing unit will 
process the stored digitized information so as to produce values of 
5 temperature and/or emissivity associated with every one of the sites 32 of the 
filter 30, for further composing a descriptive map of temperatures 
(emissivities) of the corresponding segments of the object. 

Geometric shape and parameters of the spatial filter 30 may be 
selected according to each specific application. 

io Figs. 5 and 6 schematically illustrate a particular example of the 

inventive method for obtaining temperature/emissivity maps of a selected 
region, with providing of further pattern recognition based on the obtained 
data. The particular example described below focuses on the detection of 
underground structures. It should be noted, however, that it applies mutatis 

1 5 mutandis to the detection of on ground as well as to above ground targets. 

Typically, the system used for passive remote sensing of the 
Earth's surface is a multi-spectral radiometer designed to operate in some 
convenient atmospheric window, e.g., between 8 and 12 p. The configura- 
tions and operation of such systems is well known to those skilled in the art 

20 (see. e.g., various articles appearing in: Proceedings of the First International 
Airborne Remote Sensing Conference and Exhibition, 12-15 September, 
1994, Strasbourg, France) and, therefore, will not be expounded upon herein. 
Similarly, the operation of pyranometers, for measuring optical radiation from 
the sky, is well known to those skilled in the art (Kondratyev. 1969. Radiation 

25 in the Atmosphere, 62, Academic Press) and, will also not be expounded 
upon herein. 

The detection of underground structures is performed by scanning 
the area of the Earth's surface under investigation a number of times, and 
extracting from the recorded data physical information characterizing heating 
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and cooling processes that take place not only at the Earth's surface but also 
below the Earth's surface. To this purpose an equation of thermal balance is 
used that contains information regarding heat transfer processes above and 
below the Earth's surface, through a set of therrnophysicai constants. The 
5 information obtained is then processed by suitable statistical methods in order 
to detect regions having different thermal properties, in the region selected for 
scanning. These therrnophysicai inhomogeneous regions are classified in a 
pattern recognition process in which the underground structures are identified. 



io pre-processing stage, in which temperature, emissivity and characteristic 
maps (maps of the therrnophysicai constants and statistical derivatives 
thereof) of the selected region are produced and a pattern recognition stage in 
which the underground structures are detected. 



15 scans the area of the Earth's surface under investigation at, at least one pixel at 
a time. The radiation received from each pixel by the radiometer, and by the 
pvranometer which receives radiation from the sky, serves as input for 
pre-processing software which computes characteristic maps of the selected 
region based on the Planck's equation, the albedo equation and the energy 

20 balance equation. 

A series of scans of the selected region are performed during a 
twenty- four hour period and the therrnophysicai constants are extracted from 
the energy balance equation (10): 



The detection process is divided into two main stages: a 



To effect the pre-processing stage a multi-spectral radiometer 
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-—J -T—dri (10) 

where. 7777; is the temperature of a pixel in the selected region at time rj. To 
is the temperature of the pixel at the time, to. of the first scan of the series. 
TJn) is the temperature of the atmosphere in the vicinity of the selected 
5 region at time 7 , Ifij) is the radiance of the sun and the sky in the visible and 
near infra red region, J( rf) is the radiance of the sky and the sun in the region 
8 - 1 2 u at time // . 
o is Stefan Bouzmaii's constant, 
P is the thermal inertia. 
io a is the coefficient of heat transfer, 
p is the coefficient of mass transfer, and 
O is the thermal flux for conduction. 

The energy balance equation (10) contains two functions, /and 
g that will be described below, and four unknown therrnophysical constants, 
15 a p Pand O that are to be computed at each pixel. Since, in the example 
chosen, of an underground structure, there are four unknowns, four equations 
are required to solve for them. The four equations are obtained by scanning 
the selected region five times during a twenty-four hour period. The first scan, 
perfonned at time t=0 gives rise to zero on the left hand side of equation (10) 
20 and is not used. The remaining four scans give rise to a set of four linear 
equations which can solved for the unknown therrnophysical constants. Since 
the computations are perfonned for each pixel of the selected region the result 
is a set of maps, an a map, a /?map, a P map and a Q map. 

As is well known to those skilled in the art, there are a number of 
25 possibilities for the choice of the function / for a discussion of various 
possible choices see: Jakob. 1965. Heat Transfer. John Wiley & Sons. For 
example, for low wind speeds an expression for this function is given on page 
473. Eqn. (22-91), whereas for high wind speeds an appropriate expression is 
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given on page 555 by Eqn. (26-29)). In the present example, the following 

choice for the function / is made: 

f = ( 0. 064 + 0. 032 V(r( ) os ) ( Ta (n) - T(rj) ) 

where. V( if) is the wind speed over the selected region at time t = ?] and the 
5 other variables are as defined above. This expression for the function / is 
valid for wind speeds less than 40 kilometers per hour. 

There are also a number of possibilities for the choice of the 
function g (see Rosema, 1975, "Simulation of the Thermal Behavior of Bare 
Soils for Remote Sensing Purposes" in: Report 11. Netherlands 
io Interdepartmental Working Community for the Application of Remote 
Sensing Technology, Delft, The Netherlands), in the present example the 
following choice is made: 
g('D = (n - i n /' e 

where, /„ is the time at which water starts to evaporate from the ground, at 
15 the selected region, due to heating by the sun. This time is normally defined 
as the time at which the sun makes an anele of 20° to the horizon during 
sunrise. y s and ^are constants that are determined empirically by performing 
controlled measurements for given families of ground types. Once deter- 
mined, the values of /,and /, can be stored and used as required. In general 
20 the values of y t and y 2 depend on the humidity of the ground. It has been 
found that they have one substantially constant value, for relative humidities 
less than 12%. and another substantially constant value for relative humidities 
greater or equal than 12% and less than 30%. For example, for sand-based 
earth y = 4.76 and y : = 0.72 for relative humidities of the ground less than 
25 12%. and y = 4.61 and y 2 = 0.748 for relative humidities of the ground 
greater or equal than 12%. Detection of underground structures for relative 
humidities of the ground above 30% and for rainy conditions is not possible. 
Information on rain, clouds, the atmospheric temperature and wind speed is 
obtained from meteorological stations. 
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The computation of the therm ophysical constants is summarized in Fig. 5. 
The radiation measured by the multi-spectral radiometer (block 41) along 
with relative humidity and atmospheric temperature readings (block 42) is 
conveyed to the computational process (block 43) which uses the 
5 above-described equation (4) and an albedo equation to compute the 
temperature (T), emissivity (s) and albedo (A) maps, which are stored in the 
storage means (block 44). After the fifth scan in the twenty four hour period 
has been made^ the T,£- and A maps from the last four scans along with the 
time i a and the wind speed V are conveyed to the computational process 

to (block 45). which uses the energy balance equation, to compute the maps of 
the thermophysical constants a,/3 f P and O. 

The pattern recognition stage is initiated by a process in which an 
optimal set of characteristic maps is determined. In addition to the a,j3. P 
and O maps, maps of the average, and of the variance, of a,p, Pand Q are 

1 5 computed. The average and the variance are computed by standard statistical 
techniques (see. e.g.. Tou and Gonzalez, 1974, Pattern Recognition 
Principles. Addison- Wesley Publishing Co. Inc.). In addition, texture maps 
for each of a , /?, P and O are computed. Texture maps can be computed in a 
number of ways (for a description of the concepts of image texture see. e.g., 

20 Haralick. 1982, "Image Texture Survey", in Handbook of Statistics. Krishniah 
and Kanat. editors. Vol. 2, 399), e.g., with aid of the computation of the two 
dimensional autocorrelation function according to the following equation : 

KJP) = COSf Q)fP), 

where, p is a translational distance between pixels and ^ is a translational 
25 angle between the pixels, defined as the angle between the horizontal axis and 
the line joining the origin of the coordinate system to the pixel under 
consideration. ^ and are constants of the autocorrelation function. As <j> 
varies from 0° to 360° a* and ^ take on a set of possible values, the 
parameters that describe the set of possible values of and are the 
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parameters that characterize the texture. Since the texture assists in the 
classification process it is important to have a fairly simple representation for 
the set of values it comprises. To this end, the envelope of the set of maximal 
\ allies of x$ and is chosen to be an ellipse. In order then to characterize 
5 the texture, the parameters of the ellipse of the maximal values of x$ and ^ 
have to be computed. An ellipse is described by three parameters, its minor 
axis, its major axis and the orientation of the ellipse in the plane of its 
definition. All in all. to characterize the texture, seven parameters are 
required, three that characterize the ellipse of the maximal values of x* . three 

i o that characterize the ellipse of the maximal values of co$ and one for the angle 
between the axes of the two ellipses. Hence, seven texture maps are computed 
for each of the thermophysical constants a,/3,P and 0. Hence, initially 
twenty eight texture maps are required. The maps of a,p, P and O, and of 
their average and variance, and the texture maps of a, fi, P and Q define a 

15 set of characteristic maps of the selected region. 

Areas of the selected region containing underground structures 
have surface temperatures that are very slightly different from those of the 
surroundings. These very slight temperature differences give rise to 
temperature maps of very low contrast in which it is difficult to distinguish 

20 between normal ground areas and those covering underground structures. The 
characteristic maps of the selected region serve to augment the temperature 
maps and to increase the contrast between these two types of grounds. Since 
there are a large number of maps and a large number of pixels per map, 
statistical methods are used for associating, or classifying, the pixels of the 

25 selected region with various backgrounds and targets therein. 

In the classification process the selected region is divided into 
classes. Initially the number of classes is as the number of pixels, however, in 
order for the pattern recognition process to be successful, the number of 
classes should be representative of the number of types of different 



WO 99/27336 



PCT/IL98/00568 



-32- 

backgrounds and targets, i.e., enquiry sectors, in the selected region. The 
classification process is carried out in two stages; the learning stage and the 
association stage. In the learning stage the number of classes in the selected 
region is determined, whereas in the association stage each of the pixels in the 
5 maps is associated with a class. To perform the classification a suitable 
statistical function describing the distribution of pixel values is chosen. For 
example, but by no means binding, a sum of weighted Gaussians, i.e. a 
poly-Gauss, can be chosen: 
P(s ) = Y&> NifSi.s ) 
iO where: 

Ni( Si , s ) is the i lh Gaussian in the sum, 

Wi is the weight of the i lh Gaussian in the sum, 

s is a vector whose components are the elements of the set of characteristic 

maps of the selected region, 
1 5 m is the number of classes, and 

Si represents the parameters of the i Ul Gaussian in the sum (for the properties 

of Gaussian distributions see - Tou and Gonzalez, ibid.). 

The learning stage consists of determining the values of m, co,, 

and ( >v or. if a different distribution is chosen, the equivalent parameters, co, 
20 and Si are chosen such that the poly-Gauss distribution density is maximal for 

a given vector s . The value of m is determined by including in the weighted 

sum only those terms for which : 

co/Acoi< v. 

where, v is chosen optimally by trial and error to include enough terms in the 
25 sum so that the pattern recognition process is successful, whilst on the other 
hand it should be as small as possible so as to keep the computational 
complexity to a minimum. The error in computing co,, A<y,-, is determined 
using Fisher's matrix (Rao, 1973. Linear Statistical Inference and its 
Application. Wiley). The value v=2 has been found to be optimal in a 
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number of pattern recognition experiments that where conducted. 

The association stage is performed by computing for each pixel, 
say the i'\ the probability. p ik , that it belongs to class, k, for each class in the 
selected region, using, e.g., Bayes' equation (see, Rolston, Principles of 
5 Artificial Intelligence and Expert Systems Development, McGraw-Hill. 99). 
The decision for associating the i th pixel with the k th class is made by 
determining the maximal value of 

P:A = Pa-- / A P ik , 

for a fixed value of i as k runs over all the classes. The error, Ap iki in 
io computing the probability. p ik , can be determined using Fisher's matrix (Rao. 
1973. ). The constant y, is determined empirically in controlled experiments 
by finding the minimum set of characteristic maps that enable target detection 
to a required degree of accuracy. For y = 0 the full set of characteristic maps 
has to be included in the pattern recognition process giving rise to maximum 
15 complexity in the computational process. In a certain set of controlled 
experiments it was found that J <y <3 . Invariably, the choice y-2 has been 
found to give good results for a minimal computational load. 
In the manner described, the class to which each pixel belongs is determined, 
hence effecting the association stage and thereby completing the classification 
20 process. 

Having determined the value of /, a quality factor, 0(s is 
defined in order to effect the pattern recognition process: 
Q(s )=P- yAP 

where. P is the average probability of associating each element of a set to its 
25 class, and can be computed using Bayes theorem. AP . is the error in 
computing the probability P , and can be computed using Fisher's matrix. The 
constant y. is the same as that discussed in the foregoing. The pattern 
recognition process entails determining that set of characteristic maps that 
makes O (s ) maximal. In other words the value of Q (s ) is a criterion for 
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the quality of the classification of the selected region. The process then is 
carried out as follows; choose a set of characteristic maps, perform the 
classification process and compute Q(s ). If 0(s ) is not close to the value 
L then choose another set of characteristic maps, perform the classification 
5 process again and re-compute 0(s ). If the value of 0(s ) is close to L 
then an optimal set of characteristic maps has been found and the classifica- 
tion process is declared successful. If the value of Q(s ) is not close to 1. 
then the process is repeated until the classification process can be declared 
successful. An underground structure is then an anomaly in the selected 

io region, i.e.. any area in the selected region that does not belong to a defined 
class, is then declared to be an underground structure. Normally, the partem 
recognition process starts off with a very large set of characteristic maps and 
the computational load is very heavy. In order to partially reduce the 
computational load, the pattern recognition process is performed with a 

15 priorly optimized set of characteristic maps. This optimization is achieved by 
performing a least squares calculation (see - Lawson and Hanson, Solving 
Least Squares Problems. Prentice-Hall) on Q(s ) using a weighted set of 
characteristic maps in order to find the relative weights of the characteristic 
maps in their contribution to Q (s ). Those characteristic maps having a small 

20 weight are discarded. 

Sometimes there may be prior information as to the form of the 
underground structures. In such an event the classification process can be 
enhanced and the anomalies in the selected region can be sorted by 
geometrical form, making the identification of an underground structure much 

25 more reliable. 

The pattern recognition process is summarized in Fig. 6. The 
complete set of characteristic maps is initially optimized at stage 46 using 
least squares. The initially optimized set of characteristic maps is then used to 
classify the pixels into classes (block 47), by performing the learning and 
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association stages as described above. The quality factor is then computed, at 
stage 48. and its value is checked, stage 49. If the value of the quality factor 
is not maximal then a new set of characteristic maps is chosen from the 
initially optimized set of characteristic maps (process 50). Processes 47, 48, 
5 49 and 50 are continually performed until an optimal value for the quality 
factor is obtained, at the stage pattern recognition is achieved (block 51). The 
pattern recognition process can also include an additional geometrical 
classification process for determining the shape of the detected object. This is 
particularly useful when detecting underground objects of a known shape. In 

10 general, the anomalous region defining the underground object may be larger 
than the object to be detected, in which case distribution of the anomalous 
pixels can be studied for the known geometrical trends (for standard methods 
see Tou and Gonzalez, ibid.). 

Fig. 7 schematically illustrates a system for obtaining 

15 temperature/emissivity maps of a selected object and further recognition of 
segments thereof based on the obtained maps and some additional parameters 
(such as thennophysical constants and metheorological data). An object 60 
(for example, a selected region on or above the earth) emits electromagnetic 
radiation on to an optics 62 via which it is transmitted to a multi-channel 

20 sequential filter 64 (for example, a rotational filter). The filter 64 produces N 
radiation portions in N respective spectral bands in sequence. Intensities of 
the obtained N radiation portions are then detected by a detector 66. "N" 
analogous readings of the electromagnetic radiation obtained by the detector 
66 are further converted into the digital form and transmitted to a processing 

25 unit 68 where they are stored for further processing. M*N obtained digitized 
readings of the radiation emitted from the object 60 are processed in the unit 
68 as a system of equations built according to formula (1). Values of the 
temperature (and/or the emissivity parameters or values), being solutions of 
these equations, are further transmitted to a pattern recognition processor 72 
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for building respective temperature and emissivity descriptive maps, per 
pixel, i.e. per segment of the object. However, data on the temperature and 
emissivity (if applicable) can be directly transmitted to display means 78 ; 
for optional indication. The pattern recognition processor 72 also receives 

5 meteorological and thermophysical data schematically shown as block 80, 
and processes the temperature and emissivity data received from the 
processing unit 68 to produce a file 74 being a contrast descriptive map of 
the object and a file 76 being a classified image thereof. The contrast 
descriptive map 74 is a characteristic map based on a combination of the 

10 temperature and emissivity descriptive maps, and the classified image 76 is 
a result of the pattern recognition processing based on the map 74 and the 
data 80. The display means 78 serves for graphical representation of the 
obtained files 74 and 76. The system illustrated in the drawing is adapted to 
obtain information on the electromagnetic radiation from the object 60 with 

15 the aid of a movable platform (such as a helicopter), and to process it. To 
this end. the control unit 70 of the processing unit 68 may receive data 
(schematically shown as block 82) on the velocity of the platform's 
movement with respect to the object. The unit 68 is provided with software 
means capable of registering the descriptive maps of temperature and 

20 emissivity with the platform movement. All the above-described operations 
are synchronized by a control unit 70 of the processing unit 68, which 
generates M successive synchronizing signals per each segment of the 
object 60. Each of these M signals simultaneously activates optics 62. filter 
64. detector 66 and processing unit 68 for performing and recording each of 

25 the M measurements of the radiation in N spectral bands in sequence. 
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CLAIMS: 

1 . A method of determining temperature of an object by remote sensing, 
comprising: 

- acquiring electromagnetic radiation data from the object in N spectral 
5 bands, 

- deriving temperature value from said data using a system of equations, 
each describing the electromagnetic radiation that is emitted from the object 
in said band as a function of temperature and emissivity of the object: 

the method being characterized in that, 
io - said electromagnetic radiation is measured M times, each time in N 

said spectral bands, thus obtaining N*M readings, 

- a system of N*M equations is formed, each describing the 
electromagnetic radiation that is emitted from the object in said band as a 
function of temperature T and emissivity £ of the object, wherein the 

15 emissivity e is a function of at least two arguments being temperature T and 
wavelength X: said system comprising N*M readings of the electromagnetic 
radiation and M+P unknowns which include M values of the object's 
temperature corresponding to said respective M measurements of the 
electromagnetic radiation, and P parameters, wherein P > 2: said parameters 

20 belonging to a selected mathematical form of the emissivity function; 

solving said system of N*M equations, thus deriving temperature 
of the object at any one of the M moments of said electromagnetic radiation 
measurement, wherein M is a natural number of the set [1. 2. 3...]. 

25 2. A method of determining emissivity function of an object by remote 
sensing, comprising 

- acquiring electromagnetic radiation data from the object in N spectral 
bands by measuring it M times, each time in N said spectral bands, thus 
obtaining N*M readings. 
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- forming a system of N*M equations, each describing the electromag- 
netic radiation that is emitted from the object in one of said bands as a 
function of temperature T and emissivity e of the object, wherein the 
emissivity £ is a function of at least two arguments being temperature T and 

5 wavelength X: said system comprising N*M readings of the electromagnetic 
radiation and M+P unknowns including M values of the object's temperature 
corresponding to said respective M measurements of the electromagnetic 
radiation, and P parameters, wherein P > 2; said parameters belonging to a 
selected mathematical form of the emissivity function: 
io - solving said system of N*M equations to derive said P parameters, 

thereby restoring the emissivity function of the object. 

3. The method according to Claim 1 or 2, comprising: 

- selecting said emissivity function to take constant values within 
5 respective spectral bands thus considering P equal to N, and 

- forming said system of N*M equations to comprise N*M readings of 
the electromagnetic radiation and M+N unknowns; said M+N 
unknowns including M values of the object's temperature corre- 
sponding to said respective M measurements of the electromagnetic 

o radiation and N values of emissivity for said respective N spectral 

bands. 

4. The method according to any one of the preceding claims, wherein each 
reading obtained from said M measurements of the electromagnetic radiation 

5 in N spectral bands is digitized and stored in a computer memory for farther 
computerized processing, said processing comprising substituting the 
digitized M*N readings of the electromagnetic radiation into said M*N 
equations, and solving thereof. 
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5. The method according to any one of the preceding claims, wherein said 
measurements are provided in the infrared spectrum. 

6. The method according to any one of the preceding claims, wherein said 
5 equations are integral. 

7. The method according to Claim 6, wherein said equations are Fredholm 
equations relating to emissivity. 

10 8. The method according to Claim 6 or 7. wherein each of said M*N 
equations describes the electromagnetic radiation that is emitted from the 
object in one of said bands as a function of the temperature and the respective 
emissivity of the object, the emissivity, in turn, being a known function of 
temperature, wavelength and a known number of unknown parameters: said 

15 parameters being introduced as additional unknowns into said system of 
equations for further obtaining their values upon solving the system. 

9. The method according to Claim 8. wherein said known function is 
represented as a sum of a known number of terms of two other functions, a 

20 first of which being dependent only on temperature and a known number of 
unknown parameters, and a second being dependent only on wavelength. 

10. The method according to Claim 8, wherein said known function is 
represented by a polynomial expansion having K unknown coefficients, said 

25 coefficients being introduced into the system of equations as unknown 
parameters to be determined among said unknowns upon solving said system. 

1 1 . The method according to any one of claims 1 . 8 to 10. for non-contact 
temperature measurement of heated semiconductor materials and samples 
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during manufacture thereof. 

12. The method according to any one of the preceding claims, intended for 
deriving temperature and/or emissivity descriptive maps of the object, 

5 wherein the temperature and/or emissivity is determined with respect to a 
pre- selected plurality of representative segments of the object, and mapping 
each of the determined values to a respective pixel of said descriptive map. 

13. The method according to Claim 12, wherein the object comprises a 
piuraiity of sub-objects having different temperatures and/or emissivities and 

to dislocated in a selected region. 

14. The method according to Claim 13, for detection said sub-objects with a 
background of the ground and surrounding atmosphere, comprising the steps 
of: 

1 5 - acquiring meteorological data indicative of climatic conditions of the 

selected region on the ground and/or the surrounding atmosphere: 

- recording the acquired electromagnetic radiation and meteorological 
data by means of suitable data recording means; 

- storing the so-recorded data in suitable data storage means and 
20 deriving from the so-stored data, descriptive maps of the selected 

region: and 

- feeding data of said descriptive maps into pattern recognition 
processor means to classify the pixels of said maps of said selected 
region, whereby the desired object detection is achieved. 

25 

15. The method according to Claim 14. wherein said meteorological data 
are obtained from meteorological stations. 

16. The method according to Claim 14 or 15 ; comprising providing a single 
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scanning pass of the selected region (i.e. M=I) for the detection if a 
pronounced temperature and/or emissivity differential is revealed for 
representative segments of the object. 

5 17. The method according to Claim 14 or 1 5, wherein said descriptive maps 
describe temperature, emissivity and albedo, wherefrom a set of characteristic 
maps of the selected region are derived by incorporating suitable reference 

data. 

io IS. The method according to Claim 17, wherein an optimal set of 
characteristic maps of the selected region is chosen. 

19. The method according to any one of Claims 14 to 18. wherein the pixel 
classification comprises defining a number of classes and associating the 

1 5 pixels in the selected region with said classes. 

20. The method according to any one of Claims 14 to 19, wherein a quality 
function is defined and the pattern recognition process entails determining 
that set of characteristic maps that makes the quality function maximal. 

20 

21. The method of any one of Claims 14 to 20, wherein the acquisition of 
electromagnetic radiation is performed by near range scanning. 

22. The method of any one of Claims 14 to 20, wherein the acquisition of 
25 electromagnetic radiation is performed by remote scanning. 

23. The method of any one of Claims 14 to 22. wherein the scanning of the 
selected region is performed with at least five ground oriented sensors 
operating, respectively, in the 3-4 u , 4-5 ^ 8-10 u , 10-12 u. and the near 
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infrared and visible spectral bands, and there is provided at least one further 
sensor for measuring electromagnetic radiation from the space above, 
operating in the 0.4-2 u spectral band. 

5 24. The method according to any one of Claims 14 to 23, wherein said 
classifying processor means in which the pixels of the characteristic maps of 
the selected region are classified into classes, makes use of a poly-Gauss 
expansion to describe the statistical distribution of the pixels and Bayes' 
method to sort the pixels into classes. 

10 

25. The method according to Claims 14 to 24, wherein said pattern 
recognition processor means also include additional geometrical classification 
means for determining the shape of the detected object. 

15 26. The method according to any one of Claims 14 to 25. adapted for real 
time temperature and/or emissivity mapping of a plurality of selected regions 
to be observed, wherein the step of acquisition of said electromagnetic 
radiation data is provided in real time by equipment mounted on a movable 
platform while displacing along said plurality of selected regions: and 

20 wherein the steps of recording and storing the acquired electromagnetic 
radiation and meteorological data and further deriving, from the so- stored 
data, descriptive maps of the selected region are registered with the platform 
displacement. 

25 27. The method according to Claim 26. comprising automatic recognition of 
objects positioned along or on the course of the moving platform. 

28. The method according to Claim 27, additionally comprising initiating an 
alami and automatic adjusting the course of the vehicle whenever an obstacle 
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is recognized. 

29. A system for remotely determining temperature and/or emissivity 
function of an object, said system comprising: 
5 scanning means for acquiring electromagnetic radiation data from 

the object M times, each time in N spectral bands, thus obtaining N*M 
readings of the electromagnetic radiation, wherein M is a natural number of 
the set [L 2, 3 ...]; 

means for recording and storage said N*M readings in the digital 

10 form: 

computational means for processing said N*M readings of the 
electromagnetic radiation as a system of N*M equations comprising M 
unknown values of temperature and at least P unknown parameters, wherein 
P > 2; said parameters belonging to a selected mathematical form of the 
is object's emissivity function (e) of at least two arguments being temperature T 
and wavelength \\ 

display means for displaying one or more of said unknowns 
obtained upon solving said system of equations, or derivatives therefrom. 

20 30. The system according to Claim 29, wherein said scanning means are 
capable of acquiring the electromagnetic radiation from a plural ity of 
representative segments of the object, and said computational means are 
capable of processing the acquired data, for deriving and displaying 
temperature and/or emissivity descriptive maps of the object. 

25 

3 1 . The system according to Claim 29 or 30. wherein said scanning means 
comprise a plurality of sensors capable of simultaneously measuring electro- 
magnetic radiation in a number of spectral bands. 
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32. The system according to Claim 29 to 30, wherein the scanning means 
comprise a spatial optical filter divided into a number of sites for 
simultaneous acquiring the electromagnetic radiation from said plurality of 
respectively located segments of the object; each said site having N zones 

5 responsible for splitting the radiation acquired by the site into N spectral 
bands, thereby acquiring the electromagnetic radiation in all N spectral bands 
simultaneously. 

33. The system according to Claim 29 or 30, wherein the scanning means 
comprise a rotatable optical filter for acquiring electromagnetic radiation in 

io said N spectral bands in sequence. 

34. The system according to any one of the Claims 29 to 33, wherein said 
scanning means are capable of remote scanning. 

15 35. The system according to any one of the Claims 29 to 33. wherein said 
scanning means are designed for close range scanning. 

36. The system according to any on of Claims 29 to 35. wherein said 
scanning means comprise at least five ground oriented sensors operating, 

20 respectively, in the 3-4jj., 4-5 ju, 8-10 u, 10-12 |i and the near infrared and 
visible spectral bands and at least one space oriented sensor operating in the 
0.4-2 fi spectral band. 

37. The system according to any one of Claims 29 to 36. capable of object 
25 detecting in a selected region having a background including the ground 

and/or the surrounding atmosphere, further comprising: 

- means for acquiring meteorological data indicative of climatic 
conditions of the selected region and the surroundings; 

- means for recording said meteorological data; 
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- data storage means for storing the recorded meteorological data; 

-computational means for deriving from the stored 
metheorological data and the data on the electromagnetic radiation, 
descriptive maps of the selected region; and 
5 - pattern recognition processor means for classifying the pixels of 

said maps. 

38. The system according to Claim 37, wherein said means for acquiring 
meteorological data are associated with meteorological stations. 

10 

39. The system according to Claims 37 or 38, wherein said computational 
means are capable of deriving from said descriptive maps a set of 
characteristic maps. 

15 40. The system according to Claim 37 to 39, wherein said processor means 
comprise optimizing means for selecting an optimal set of characteristic maps 
which is subjected to pixel classification. 

4 1 . The system according to any one of Claims 37 to 40, wherein the pattern 
20 recognition processor means comprises pixel classification means in which a 

number of classes is defined whereby the classified pixels are statistically 
associated with the classes. 

42. The system according to any one of Claims 37 to 4 1 , wherein the pattern 
25 recognition processor means is capable of functioning based on a poly-Gauss 

expansion for determining the statistical distribution of the pixels, and means 
based on Bayes's method for sorting pixels into classes. 

43. The system according to any one of Claims 37 to 42. wherein said 
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pattern recognition processor means comprises quality factor evaluation 
means. 

44. The system according to any one of Claims 37 to 43, wherein the pattern 
5 recognition processor means includes additional geometrical classification 

means for determining the shape of the detected object. 

45. The system according to any one of Claims 37 to 44, wherein the 
scanning means and the data recording means are positioned on any movable 

io platform, whereas the data storage means, computational means, classification 
means and pattern recognition means are placed on the ground. 

46. The system according to any one of Claims 37 to 44. being portable and 
mountable on a movable platform. 

15 

47. The system according to Claim 45 or 46, wherein said computational 
means are adapted to register the recorded and stored data of the acquired 
electromagnetic radiation with the platform's movement to obtain a real time 
mapping and pattern recognition. 

20 

48. The system according to any one of Claims 45 to 47, wherein said 
pattern recognition processor means are capable of automatic recognition of 
objects positioned along or on the course of the moving platform. 

25 49. The system according to Claim 48, additionally comprising means for 
initiating an alarm. 

50. The method according to any one of Claims 14. 15. 17 to 26. for 
detection of underground structures, comprising extracting from the recorded 



WO 99/27336 



-47- 



PCT/IL98/00568 



data physical information characterizing heating and cooling processes that 
take place at the ground surface and below it by applying an equation of 
energy balance and a set of thermophysical constants, processing said 
information for detecting regions having different thermal properties, in the 
5 region selected for scanning, and classifying thus revealed thermophysical 
inhomogeneous regions by said pattern recognition processor means, thereby 
identifying the underground structures. 

Dl. The method according to Claim 50 ; comprising derivation of the 
io descriptive and characteristic maps based on Planck's equation, describing 
the electromagnetic radiation that is emitted from and reflected by the 
ground. 

52. The method according to Claim 51, wherein said characteristic maps 
15 comprise maps of thermal inertia, thermal flux, coefficient of heat transfer 

and coefficient of mass transfer. 

53. The system according to Claims 37 to 47, wherein said pattern 
recognition processor means are capable of extracting from the recorded 

20 data physical information characterizing heating and cooling processes that 
take place at the ground surface and below it by applying an equation of 
energy balance and a set of thermophysical constants, processing said 
information for detecting regions having different thermal properties, in the 
region selected for scanning, and classifying thus revealed thermophysical 

25 inhomogeneous regions, thereby enabling identifying and displaying the 
underground structures. 
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Underground Physical Non-uniformity Detecting Method. 



The invention is related with geological recognition methods and may be used in 
mineral products searching, geological investigations, etc. The invention goal is the 
improvement of underground non-uniformity detection method by means of low 
contrast in the detection of thermo-physical properties objects. 

This method may be realized in the following way: 

Twice remote Earth surface thermo surveying (mapping) must be done on the tested 
area territory. First surveying must be done after sunset, and second - before sunrise. 
Some reference zones within tested area must be selected. These selected zones must 
be characterized by presence of differences of thermo-physical properties objects. 

Night thermo contrast inversion moment must be defined by means of a few 
methodical surveyings. That is the moment, when different thermo inertia mineral 
bodys will be equal on temperature during the night. The following data processing 
must be done depending on above mentioned thermo inversion effect presence or 
absence. If thermo invertion effect present, substraction of two thermo maps must be 
used in following data analysis. If thermo inversion effect absent, adding of two maps 
must be used. If resultant thermo map has an anomalous, this map may be used for 
underground irregularity detection. 

Optimisation for remote recognition season choice must be done with night thermo 
inversion moment for different thermo inertia bodys taking into account. If thermo 
inversion effect is present, surveying must be done at the time, when the inversion 
moment is maximum removal from sunset moment. If thermo inversion effect is 
absent, surveying must be done at the time of maximum minerals temperature 
difference, discovered on reference zone during the night. 

This method provides the low contrast thermo non-uniformity detection. 
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HUH ManOKOHTpaCTHbDC no TenJTO(pH3KMeC- 
KHM CBOHCTBSM 06T>eKTOB.. GllOCOti peajTH- 

syercfl nyreM flByKpaTHoro HSMepe hmh 
TeruroBoro nonw s reveHHe oahoA- howh 
npH onTHMH3aiyiH ce30HHoro speMeim 
o n pefl ejie n h_h m HHTe pnpcTa u,hh peayjib- 
THpyiome ro no/i« b 3sbxcmmocth ot no/io- 

XeHHfl TOWKH HOMHOfl HHBfipCHH . 
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HpOH36DflCTBeH 



HO-no^rpa^H^ccKoe npeflnpM«Tne, r.y^ropofl , yn.npocKTi.aH, 




1 



1383259 



M3o6peTeHHe othochtch k reoTep- 

MHHeCJCHM MeTOflaM pa3BeflKH H MOieT 

BbiTfa HcnoJit>3 0BaHO npH noHCKax none3- 
HbK HCKonaeMbix, HHJKeHepHo-reojiorHMec- ^ 

KHX H3bJCKaHHHX H T.fl. 

UenB H^oSpeTeHHfl — pacnmpeHHe ao3— 
MoistnoCTeH cnocoGa oGHapysceHHH noflseM— 
HbK HeoAHopoflHOCTeH aa CMeT BblHBJie- 

HHH MflJlOK O H TP a C T HblX FtO TerUIOlpHSHHeC- 10 
KHM CBOftCTgaM 06l»eKTOB. 

Cnocoo peaJtHayeTCs cjieAyiomHM 06- 
Ha HCCJieAyeMOH reppHTopHH (no-riHro- 

He) . npOBOflHT AByKpaTHyiO AHCTaHIJHOH- "15 

Hyw TemioByio cteMKy noBepxHocTH 3eM- 
jih. IlepBoe H3MepeHHe ocymecTBJinioT .no- 
cjie 3axoAa, a BTopoe - b npeApaccBer- 
Hoe BpeMH nepeA bocxoaom ConHiia. 

B npeflejiax nonHTOHa BbiSHpaioT 3Ta- 20 
noHHbie ynacTKH, xapaKTepH3yiomHecH Ha— 
jth^hcm nopoA c p a 3 jih*i HbiMH Teruio4>H3H- 

HeCKHMH CBoftCTBaMH, Ha Bbl6paHHbDC yHil~ 

ctk3Xj HanpHMep, c noMombw pexHMHbix 
Ha 6 jiioa e jmi5 , . onpeAeJiawT BpeMH hovhoh 25 
HHBepcHH TeMnepaTypKbix kohtpsctob t 
T.e*. BpeMH, KorAa xeMnepaTypw nopoA 
c pasJiHKHofi TeiuioBoft HHepAHen CTaHo- 
BHTCH pa,BHfaiMH.B TeneHHe how, 

B 3 a BHC HMO C TH OT HaJW^UiH HJIH OT" . 30 
Cy.TCTBHfl HOIHOft HHBGpCHH o5pa6oTKy 

AaHHbrx, nany^ieHHbix no one npoBeAeHKH 
AsyxpaTHOH AHCTaHUHOHHoft TennosoH 
cteMKH- Ha .noJiwroHe, hpohsboaht cjie— 
AyramKM ofipa30M. 35 

IIpH HaJTHTOH HOMHOfl KHBfipCHH TeMnfi- 

paTypHbix KOHTpacTpB pe3yAtTHpyiomee 
TeiuioBoe none onpeAejiawT no pasHocTH 
nojiyHeHHbix TeMnepaTyp, a npH ee ot- 
cyTCTBKH - no hx cyMMe. no aHoMajib- 4Q 
hum 3 h a^s hhhm pe3yjii>THpyiomero Teano- 
Boro noJiH BhiHOCHT cyatfleHHe o Ha/iHwii 
noAseMHbix HeoAHOpOAHocTeH • 



Bw6op onTHMa/ibHoro ceioHHoro nepH- 

OAa AHH npOBeAeHHH A**CTaHUHOHHbiX 

H3MepeHHH Taoe npOH3BOflflT c y^eTOM 
BpeMeHH hohhoh KHB£pcHH TeMnepaTyp- 
Hbix KOHTpacToB nopoA c pa3JTHMHOH Ten— 

JIOBOH HHepAHeft. 

. IIpH HajlHHHH HOWHOH HHBepCHH .H3Me- 

peHHH npoBOflflT b Te^teHHe nepHOfla, 
Korfla BpeMH hhbspchh MaxcHMajifcHo yAa- 
neHO ot t iaxoAa CojtHua., npH otcvtctbhh 

H0MHOH HHBepCHH H3MepeHHfl OCymeCTB- 

jihiot b nepwoA HaHdoji.&oiHX nepenaflOB 
TeMnepaTyp nopOA Ha DTanoHHOM y^acT- 
Ke b TeqeHHe ho^im. 

Hcnojit30BaHHe cnoco6a o6ecnewsaeT 
BWHBJteHHe MajxpKOHTpacTHbDC no Tenno- 

(t»H3HMeCKHM CBOHCTB3M HeOAHOpOflHOCT ePl . 

G> o p m y n a h a o 6 p e t e h h h 

CnocoC oOHapysceHKH noA3eMHbtx He- 
OAHOpoAHocTeii c pa 3-nHHHMMH Tenno<$>H3H— 

leCKHMH CBOACTBSMK f - BKJIIOHaKHUHH ABy** * 

KpaiHoe ' npoBeAeHHe AHCTaH^oHHofl Ten- 

JIOBOA . C-faeMKH TIOBepXHOCTH 3eMJIH, oaho 

H3 KOTopbix npoH3BOflHTCfl b npeApac-. 
CBeTKoe BpeMH 9 o6pa6oTxy nojiyMeHHbix 
A an HbK h cy^qeHHe o nasimnn HeoAHO-- 

pOAHOCTefi, OTJl-HMaiomHHC H 

TeM, <tTo, c Aejxfcio pacumpeKHfi. bo3mo3k— 
HocTeH cnocoSa aa cieT Bbbr BJie hhh Ma— 

HOKOHTpaCTHbDC rio TenJl6<pH3H*ieCKHM 

CBoficTBaM o6ieKTOB t BTopoe " HSMepemie 
npoBOAHT nocjie saxoAa conriua, onpeAe- 
ji;iiot BpeMH HosHoft HHBepCHH TeMnepa- 

TypHbJX KOHTpaCTOB nopOfl C pa3JlHMHOH 

TennoBOft HHepAHeft, npH HajiH/iHH hom- 
hoh HKBep^tiH pe3yjn>THpywmee Terniosoe 
'no/ie onpeAeJiHroT no- pasHocTH nojiy^ieH- 
Hboc TeMnepaTyp, a npH ee oTcy tctbhh - 
no hx cyMMe . 
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